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SECTION  I 


INTRODUCTION 

It  is  well  known  that  cracks  whicn  occur  at  fastener  holes  or  attach¬ 
ment  lugs  are  common  sources  of  component  failure  in  aircraft  structures. 
Reference  1  presents  a  1971  review  of  significant  USAF  aircraft  structural  fail¬ 
ures  which  revealed  that  bolt  and  rivet  holes  comprised  over  one  third  of 
the  failure  origins  considered.  In  both  the  lug  and  fastener  hole  geome¬ 
tries,  high  stresses  are  generated  next  to  the  edge  of  the  hole  by  load 
transfer  through  the  fastener  or  pin.  These  hign  stress  levels  can  quickly 
initiate  cracks  at  local  imperfections  along  the  bore  of  the  hole.  Sub¬ 
sequent  loading  may  then  cause  further  flaw  growth  and  ultimate  failure 
due  to  the  combined  action  of  fatigue  and  corrosion. 

The  objective  of  this  report  is  to  summarize  results  of  an  effort 
directed  at  characterizing  the  growth  of  fatigue  cracks  which  initiate 
midway  along  the  bore  of  a  hole  in  a  large  plate  or  attachment  lug.  The 
specimen  geometries  and  crack  configurations  of  interest  are  shown  in 
Figures  1  and  2.  The  plate  specimens  contained  open  holes  and  were  loaded 
in  remote  tension,  while  the  lug  specimens  were  loaded  through  a  pin 
placed  in  the  hole.  In  all  cases  the  crack  plane  was  located  perpendicu¬ 
lar  to  the  loading  axis. 

In  both  specimen  geometries,  semielliptical  surface  cracks,  with 
major  and  semiminor  axes  2a  and  c,  were  introduced  midway  along  the  hole 
bore  as  shown  in  Figure  2a.  These  cracks  were  extended  by  fatigue  and 
monitored  as  they  penetrated  the  specimen  surfaces  and  transitioned  into 


Uniform  Through-The-Thickness  Cracks. 


uniform  through-the-thickness  flaws.  Immediately  after  the  semielliptical 
crack  penetrated  the  left  and  right  specimen  surfaces  (see  Figure  2b),  the 
flaw  shapes  were  characterized  by  the  surface  dimensions  c^  and  cR  and  by 
the  midpoint  length  c.  Once  a  uniform  through-the-thickness  geometry  was 
reached,  the  five  measurements  shown  in  Figure  2c  were  used  to  determine 
an  average  through-the-thickness  crack  length. 

The  test  matrix  is  summarized  in  Table  1  along  with  nominal  specimen 
dimensions.  Note  that  three  plate  specimens  were  tested  with  single  embed¬ 
ded  surface  cracks,  while  one  plate  contained  two  nonsymmetric  corner 
cracks  located  at  opposite  sides  of  the  hole.  Nine  pinloaded  lug  specimens 
were  tested  with  single  embedded  surface  cracks.  Five  of  these  specimens 
were  continued  as  uniform  through-the-thickness  flaw  tests.  The  lug  speci¬ 
mens  included  four  different  widths  (W/D  =  3.0,  2.5,  2.0,  and  1.5).  Here 
W  is  the  specimen  width  and  D  is  the  hole  diameter.  Three  four-point  bend 
edge-crack  specimens  were  also  tested  to  further  characterize  the  baseline 
fatigue  crack  growth  properties  of  the  specimen  material. 


Since  direct  measurement  of  the  internal  dimensions  of  the  semiellip¬ 
tical  and  through-the-thickness  cracks  would  be  impossible  in  metal  speci¬ 
mens,  the  test  pieces  were  manufactured  from  polymethylmethacrylate  (PMMA), 
a  transparent  polymer.  The  transparent  nature  of  PMMA  allows  internal 
crack  dimensions  to  be  monitored  during  the  test,  and  provides  a  simple, 
but  accurate  means  for  measuring  both  crack  shape  and  size  changes  as  a 
function  of  applied  loading  cycles. 

Being  a  fairly  brittle  material,  flaw  growth  in  PMMA  may  be  described 
by  many  of  the  same  linear  elastic  fracture  mechanics  techniques  employed 
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for  conventional  structural  alloys.  In  particular,  both  fracture  and 
fatigue  crack  growth  rates  can  be  characterized  in  terms  of  the  stress 
intensity  factor  K.  The  following  section  describes  the  experimental  pro¬ 
cedure  and  material  characterization  tests,  while  the  experimental  results 
and  corresponding  fracture  mechanics  analyses  are  described  in  later  sec- 
ti ons . 
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Test  Matrix 


Two  nonsymmetric  corner  cracks  located  at  opposite  sides  of  hole. 


SECTION  II 


EXPERIMENTAL  PROCEDURE 

The  transparent  property  of  PMMA  was  used  in  the  present  experiments 
to  make  direct  measurements  of  the  embedded  surface  cracks.  The  methods 
used  to  characterize  the  test  material  and  to  measure  the  cracks  in  the 
PMMA  plate  and  lug  specimens  are  described  in  this  section.  Additional 
details  may  be  found  in  Reference  2. 

1 .  Material  Characterization 

All  specimens  were  cut  from  a  single  sheet  of  PMMA  and  were  located 
so  that  the  crack  plane  had  the  same  orientation  with  respect  to  sheet 
reference  axes.  The  PMMA  sheet  had  supplied  specimens  for  previous  pro¬ 
grams  (3-6)  and  was  fairly  well  characterized.  The  0.2  per  cent  tensile 

yield  strength  was  previously  found  to  be  7000  psi  and  the  fracture  tough- 

1  /2 

ness  KjC  was  990  ±  60  psi-in 

Baseline  fatigue  crack  growth  rate  tests  were  conducted  with  edge- 
cracked  beams  loaded  in  four-point  bending  as  shown  in  Figure  3.  The  cy¬ 
clic  load  was  applied  at  a  frequency  of  2HZ  with  0  <  R  £  0.05  (R  =  minimum/ 
maximum  load).  The  crack  plane  was  illuminated  with  a  fibre  optics  light 
source  and  photographed  through  the  transparent  end  of  the  specimen  at 
periodic  cyclic  intervals.  A  Nikon  F3  camera,  135  mm  lens,  and  bellows 
were  used  to  photograph  the  cracks.  A  typical  baseline  test  photograph 
is  shown  in  Figure  4. 

The  crack  photographs  were  measured  by  projecting  the  35  iron  black 
and  white  negatives  onto  a  screen.  Measuring  the  projected  images  to  the 
nearest  millimeter  resulted  in  an  accuracy  of  0.001  in.  for  the  actual 
crack  lengths.  Since  the  crack  front  is  curved,  a  five  point  average  was 
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Figure  4  Photograph  of  a  Through-The-Thickness  Edge  Crack 
in  a  Four-Point  Bend  Baseline  Specimen 
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used  as  the  crack  length  (see  Figure  2c).  ASTM  guidelines  [7]  recommend 
that  cyclic  crack  growth  increments  be  on  the  order  of  0.01  in.  for  this 
level  of  precision.  This  criterion  was  used  as  much  as  possible  when 
choosing  cyclic  intervals  between  photographs.  The  crack  length  data  for 
the  tests  are  shown  graphically  in  Figure  5.  Fatigue  crack  growth  rates 
(da/dN)  were  determined  by  a  five-point  polynomial  method  [8]  for  differ¬ 
entiating  the  experimental  data. 

The  baseline  da/dN  versus  cyclic  stress  intensity  factor  data  are 
compared  in  Figure  6  with  a  band  representing  data  reported  in  Reference  4 
for  the  same  sheet  of  material.  (The  Reference  4  data  were  obtained 
at  a  cyclic  frequency  of  1  Hz.)  Note  that  the  crack  growth  rate  is  described 
quite  well  over  the  range  shown  by  the  simple  power  law 

^=CAKm  (1) 

Here  C  and  m  are  the  empirical  constants  given  on  Figure  6.  Stress 
intensity  factors  were  computed  for  the  four-point  bend  geometry  by  the 
appropriate  expression  given  in  Reference  9. 

2.  Plate  and  Lug  Specimen  Preparation 

Steel  grips  were  bonded  and  bolted  to  the  plate  and  lug  specimens  as 
shown  in  Figures  7  and  8.  A  mirrorplaced  over  the  specimen  ends  allowed 
direct  observation  of  the  crack  plane.  Preparation  of  the  plate  and  lug 
specimens  included  five  steps:  machining,  annealing,  surface  preparation, 
bonding,  and  crack  initiation, 
a.  Machining 

The  plate  and  lug  specimens  described  in  Figure  1  were  machined 
to  proper  dimensions  from  the  same  0.7  in.  thick  PMMA  sheet  used  for 
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Figure  6 


)N  =  CCAtffJ  C  =  1.0701783*  10s1 
ms  9.21 4 13  5  6 


its  of  C  are  for  da/dN  in  IN. /CYCLE 
i  AK  in  PSI /IN. 


Baseline  da/dN  vs.  aK  data  for  PMMA 
Test  Material 
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Figure  8  Schematic  View  of  Pin-Loaded  Lug 

Specimen  Fastened  to  the  Test  Grips 
(All  Dimensions  in  Inches). 
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the  baseline  specimens.  As  mentioned  earlier,  all  of  the  specimens 
were  cut  so  that  the  cracks  would  grow  in  the  same  direction  relative 
to  the  sheet. 

b.  Annealing 

Each  specimen  was  annealed  in  an  oven  at  101 °C  (214°F)  for  24 
hours  to  relieve  any  residual  stresses  that  may  have  been  present 
after  machining.  The  specimens  were  cooled  to  50°C  (122°F)  at  a  rate 
of  2°C  per  hour  (3.6°F  per  hour).  Then  the  heat  was  turned  off  and 
the  specimens  were  allowed  to  cool  to  room  temperature  (21 °C  =  70°F) 
with  the  oven  door  closed.  Total  cooling  time  was  approximately  40 
hours. 

c.  Surface  Preparation 

The  PMMA  material  was  covered  with  protective  paper  that  needed 
to  be  removed  prior  to  bonding  to  the  metal  grips.  Since  the  material 
used  here  was  several  years  old,  the  paper  was  difficult  to  remove. 

A  flat  wood  file  proved  to  be  the  most  efficient  way  to  remove  the 
paper.  The  bonding  areas  were  then  sanded  with  120  weight  sandpaper 
for  good  bond  adhesion  between  the  specimens  and  loading  grips. 

The  viewing  end  of  each  specimen  serves  as  a  lens  and  must  be 
clear  and  flat  in  order  to  view  the  cracks  without  distortion.  The 


end  was  sandwiched  flush  between  two  flat  pieces  of  scrap  material 
and  clamped  in  place.  A  sanding  block  was  used  to  sand  the  surface 
to  a  uniform  texture  with  the  following  order  of  sandpaper:  120,  220, 
440,  and  600  weight.  The  final  step  was  to  polish  the  end  to  trans¬ 
parency  with  one  micron  aluminum  oxide  powder  on  a  moist  cloth. 
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d.  Bondi ng 

All  bonding  surfaces  were  wiped  clean  with  alcohol.  Thermoset 
103  adhesive  was  used  with  a  resin/hardener  ratio  of  2/1  which  pro¬ 
vides  a  tensile  shear  strength  of  3,100  psi  [10].  The  grip/specimen 
alignment  was  provided  by  a  gluing  rig.  Bonding  pressure  was  applied 
by  a  25  pound  lead  weight,  and  silica  grease  was  used  to  prevent  bond 
adhesion  at  the  viewing  surface  or  grip/gluing  rig  interface.  Using 
holes  in  the  steel  grips  as  guides,  holes  were  drilled  through  the 
PMMA  specimens  and  used  to  bolt  the  grips  to  the  test  specimens  (in 
addition  to  the  bonding). 

e.  Crack  Initiation  Methods 

Before  any  of  the  specimens  could  be  tested,  an  initial  embedded 
flaw  had  to  be  introduced.  The  plate  specimens  had  a  hole  diameter 
of  0.75  in.  which  gave  little  working  room  for  starting  a  crack.  A 
pointed  soldering  iron  was  used  to  make  an  indentation  in  the  center 
of  the  hole  normal  to  the  loading  direction.  This  s  :ress  concentration 
caused  a  crack  to  start  at  that  dc,i  •">.  ThtJ  tie  diameter  of  the  lug 
specimens  was  large  enough  to  allow  the  use  of  a  X-acto  knife  for 
crack  initiation.  The  knife  was  positioned  in  the  hole  and  impacted 
with  a  small  hammer  to  cause  a  precrack  to  propagate  in  front  of  the 
blade.  Late  in  the  testing,  it  was  found  that  this  method  occasion¬ 
ally  caused  some  crack  retardation  due  to  crack  tip  blunting,  and  in 
some  cases  caused  cracks  to  grow  out  of  plane.  A  hot  X-acto  knife 
was  used  to  initiate  cracks  in  lug  tests  PT5  and  PT11 . 

All  of  the  precracks  were  further  extended  with  cyclic  loads 
approximately  40-50%  larger  than  the  final  test  load.  A  frequency  of 
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3  Hertz  was  used  for  this  precracking.  Although  PMMA  is  brittle, 
and  shows  little  fatigue  crack  retardation  [3],  care  was  taken  to 
reduce  the  loading  slowly  during  precracking  to  avoid  any  retardation 
at  the  final  test  loads.  The  initial  crack  lengths  along  the  hole 
were  precracked  to  approximately  50%  of  the  specimen  thickness.  All 
of  the  final  crack  growth  testing  was  conducted  at  the  same  cyclic 
frequency  of  2  hertz  used  for  the  baseline  testing  in  order  to  avoid 
possible  rate  effects  in  the  polymer  test  material.  The  same  R  ratio 
(0  <  R  _<  0.05)  was  used  for  these  tests  as  for  the  baseline  tests. 

It  should  be  noted  that  an  attempt  at  making  small  initial  flaws 
was  made  with  the  lug  specimens.  This  was  accomplished  by  drilling 
a  smaller  two  inch  diameter  hole  in  the  lugs.  The  lugs  were  then 
precracked,  and  the  hole  was  bored  to  its  final  size  (2.25  inch) 
leaving  a  small  initial  flaw.  Since  the  small  flaws  (<0.5  T)  did  not 
yield  fast  enough  crack  growth  rates  at  the  test  loads,  and  the  pro¬ 
cedure  was  fairly  time  consuming,  this  method  was  abandoned.  However, 
it  was  used  on  lug  tests  PT4  and  PT6. 

3.  Fatigue  Crack  Growth  Testing 

A  10  kip  programnable  MTS  machine  was  used  to  load  the  specimens. 

The  grips  were  pinned  to  the  load  cell  and  the  hydraulic  piston  through 
eyebolts.  A  movable  arm  was  designed  and  used  as  a  mount  for  the  camera 
and  bellows.  The  cracks  were  viewed  through  a  front  surface  mirror  set 
at  approximately  45  degrees  in  the  upper  grips.  A  double  gooseneck  fiber 
optics  (cold)  light  source  was  used  to  illuminate  the  cracks.  The  cold 
light  source  avoided  heating  the  test  specimens.  The  testing  arrangement 
is  illustrated  in  Figure  9  for  one  of  the  plate  specimens. 
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Figure  9a  Photograph  of  Testing  Apparatus  Showing  Place  Specimen 
Mounted  in  Fatigue  Machine,  Viewino  Mirror,  Camera, 
Light  Source,  and  Strain  Indicator. 
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Figure  %  Closeup  View  of  Test  Apparatus  Showing  Camera  and  Viewing  Mirror  for  Plate  Specimen 


Although  the  specimens  were  pinned  at  both  ends,  out  of  plane  bending 
may  have  existed  if  the  specimens  were  not  loaded  through  their  axis  of 
symmetry.  All  of  the  plate  specimens  were  fitted  with  strain  gages  on 
each  side  at  the  crack  plane.  These  gages  monitorea  strain  differences 
through  the  thickness  of  the  material.  Washers  were  used  to  align  the 
specimens  to  minimize  bending.  Strain  differences  greater  than  five  per¬ 
cent  could  be  detected  visually.  Two  early  practice  lug  tests  were  also 
mounted  with  strain  gages,  and  showed  similar  results.  Lug  tests  PT4, 

PT6,  and  PT8  were  positioned  visually.  All  of  the  other  lug  tests  were 
aligned  by  using  a  dial  indicator  to  sense  horizontal  motion  below  the 
crack  plane. 

4.  Crack  Measurements 

Several  factors  are  involved  in  viewing  the  cracks  in  the  plate  and 
lug  specimens.  First,  the  viewing  surface  must  be  clear  and  flat  since  it 
is  a  lens.  Also,  the  light  has  to  shine  nearly  parallel  to  the  camera  so 
that  the  crack  reflection  may  be  observed.  One  of  the  biggest  problems 
that  was  encountered  was  the  effect  of  specimen  length  on  crack  observation. 
It  was  very  difficult  to  see  both  sides  of  the  lug  specimens  if  the  dis¬ 
tance  from  the  lens  to  the  crack  plane  was  greater  than  8  inches.  The 
shortest  specimens  that  still  gave  uniform  stress  at  the  crack  plane  were 
determined  to  be  those  with  an  unbonded  length/width  (L/W)  of  unity.  This 
was  determined  from  stress  intensity  factor  solutions  for  similar  geome¬ 
tries  that  gave  "infinite  length"  stress  intensity  factors  [9].  There¬ 
fore,  all  of  the  lug  specimens  had  unbonded  L/W's  of  one.  The  unbonded 
plate  lengths  were  11  inches.  The  bond  length  in  all  cases  was  1.5  inches, 
so  that  the  total  length  to  the  crack  plane  war  7  inches  for  the  plate 
specimens . 
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Crack  growth  was  again  recorded  by  time  lapse  photography  for  the 
plate  and  lug  tests.  A  typical  set  of  photographs  for  one  of  the  plate 
tests  is  shown  in  Figure  10.  Here  the  flaw  plane  is  oriented  so  that  the 
hole  bore  is  vertical  and  the  crack  is  growing  toward  the  right  from  the 
right  edge  of  the  hole.  Note  that  the  transparent  specimen  acts  as  a 
prism  and  the  crack  image  is  reflected  at  the  specimen  walls.  The  initial 
surface  crack  grows  along  the  hole  bore  and  Denetrates  one  of  the  plate 
surfaces  at  approximately  46,300  cycles.  Penetration  through  the  opposite 
surface  follows  shortly  after,  and  the  crack  continues  to  grow  into  a 
through-the-thickness  geometry.  (Detailed  measurements  are  reported  later.) 

The  crack  photographs  were  measured  as  before.  The  negatives  were 
projected  to  a  scale  so  that  0.001  inch  in  actual  size  could  be  resolved. 
Scaling  of  the  photographs  was  based  on  known  reference  dimensions  in  two 
perpendicular  directions.  The  hole  diameter  and  plate  thickness  provided 
convenient  reference  distances.  Slight  scribe  marks  along  the  front  and 
back  surfaces  in  the  crack  plane  were  sometimes  used  to  highlight  the 
plate  thickness  in  the  photographs. 

Since  the  embedded  flaws  eventually  transition  to  become  through-the- 
thickness  cracks,  the  crack  shapes  were  divided  into  three  categories: 
embedded,  intermediate,  and  transition  crack  lengths.  The  embedded  case 
implies  that  the  crack  is  completely  inside  the  material  along  the  bore 
of  the  hole  (see  Figure  2a).  The  intermediate  case  occurs  when  one  side 
of  the  crack  has  penetrated  the  surface  of  the  material  while  the  other 
side  is  still  embedded.  This  is  a  result  of  unsymmetrical  initial  flaws, 
material  inhomogeneity,  and  the  existence  of  some  out  of  plane  bending. 

When  both  sides  of  the  crack  have  penetrated  the  surface,  as  shown  in 


Figure  2b,  transition  is  achieved.  In  the  embedded  case,  the  hole  bore 
crack  dimension  2a,  and  the  maximum  length  of  the  crack  radiating  from 
the  hole  (dimension  c)  were  measured  as  a  function  of  elapsed  load  cycles. 
Crack  lengths  along  each  side  of  the  specimen  (c^  and  cR),  as  well  as  the 
maximum  penetration  c,  were  measured  after  transition.  In  the  intermediate 
case,  after  surface  penetration  by  one  crack  tip,  but  not  both,  the  hole 
dimension  2a  was  also  recorded.  The  crack  growth  results  are  summarized 
in  the  following  section. 


SECTION  III 


EXPERIMENTAL  RESULTS 

The  results  of  the  experiments  conducted  here  may  be  divided  into 
three  categories:  plate  tests,  lug  embedded  crack  tests,  and  lug  through- 
crack  tests.  All  of  the  plate  specimens  were  loaded  in  remote  tension, 
and  the  cracks  grew  from  the  open  holes.  There  were  four  different  lug 
sizes,  all  of  which  were  pin  loaded.  The  data  for  each  test  are  given  here 
in  graphical  format  and  are  presented  in  tabular  form  in  Reference  11.  The 
actual  specimen  dimensions  and  applied  cyclic  force  for  each  test  are 
summarized  in  Table  2. 

1 .  Plate  Test  Results 

All  of  the  actual  crack  lengths  were  calculated  by  scaling  the  photo¬ 
graph  measurements  using  known  dimensions  as  discussed  in  Section  II.  These 
crack  lengths  were  plotted  versus  elapsed  load  cycles  before  and  after  both 
sides  of  the  crack  had  penetrated  the  specimen  thickness.  Before  either 
side  of  the  crack  penetrated  the  plate  thickness,  crack  dimensions  a  and  c 
were  measured  and  plotted  (see  Figure  2).  After  one  side  of  the  crack  had 
"popped"  through  a  plate  surface,  a  new  plot  was  prepared  to  describe  the 
transition  and  growth  of  the  through-the- thickness  flaw.  This  second  plot 
presents  the  dimensions  c^,  c^,  and  c  for  the  through-crack,  and  for  refe¬ 
rence,  also  includes  the  last  5  measurements  of  the  embedded  flaw  prior  to 
free  surface  penetration. 

The  crack  growth  data  for  plate  specimens  T1 ,  T2,  T3  and  Ul  are  shown 
in  Figures  11-14.  Note  the  crack  growth  behavior  following  penetration  of 
the  specimen  thickness  by  the  embedded  cracks  in  Specimens  Tl,  T2,  and  T3 
(Figures  lib,  12b  and  13b).  Here  the  new  crack  dimensions  created  at  the 
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Table  2  Specimen  Dimensions  and  Applied  Loads 


Specimen  W  T  D  Applied  Load  P(lbs) 

Number  (inch)  (inch)  (inch)  Minimum  Maximum 


Surface  Cracked  Plate  Specimens 


T1 

7.975 

0.699 

0.750 

10 

3400 

T2 

7.970 

0.720 

0.750 

25 

3400 

T3 

7.960 

0.700 

0.750 

50 

3400 

U1 

8.000 

0.699 

0.750 

30 

3400 

Surface  Cracked  Lug  Specimens 

PT4 

6.750 

0.712 

2.250 

10 

1100 

PT5 

6.750 

0.715 

2.250 

50 

1200 

PT6 

5.625 

0.711 

2.250 

10 

1100 

PT7 

5.625 

0.714 

2.250 

10 

1000 

PT8 

4.500 

0.714 

2.250 

5 

750 

PT1 4 

4.500 

0.714 

2.250 

20 

900 

PT10 

3.375 

0.724 

2.250 

5 

500 

PT1 1 

3.375 

0.711 

2.250 

5 

500 

PT1 3 

3.375 

0.750 

2.250 

0 

800 

Through- 

-Cracked  Lug 

Specimens 

PT4 

6.750 

0.712 

2.250 

10 

1100 

PT5 

6.750 

0.715 

2.250 

50 

1000 

PT8 

4.500 

0.714 

2.250 

5 

550 

PT9 

4.500 

0.724 

2.250 

5 

550 
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Figure  11a  Fatigue  Crack  Growth  Data  for  Plate  Specimen 
T1  F<"ior  to  Crack  Transition 


FMBEDDED  CASE 


T3 
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Figure  13a  Fatigue  Crack  Growth  Data  for  Plate  Specimen 
T3  Prior  to  Crack  Transition 
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RACK  LENGTH  UN. 


EMBEDDED  CASE 


TRANSITION  CASE 


N  (CYCLES) 


Figure  14  Fatigue  Crack  Growth  Data  for  Nonsymmetric  Double 
Corner  Crack  Plate  Specimen  U1 
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plate  surfaces  (cL  and  cR)  initially  grow  very  rapidly,  but  then  slow  down 
as  a  uniform  through-the-thickness  geometry  is  reached.  These  dimensions 
then  increase  rapidly  again  prior  to  final  fracture.  This  transition  beha¬ 
vior  is  correlated  with  corresponding  stress  intensity  factor  analyses,  and 
is  discussed  in  greater  detail  in  a  separate  paper  [12], 


A  composite  plot  of  crack  aspect  ratio  a/c  versus  flaw  size  2a/T  is 
given  in  Figure  15  for  the  three  embedded  plate  tests.  Note  that  all  three 
tests  yielded  similar  crack  shapes  prior  to  penetration  through  the  plate 
surfaces.  The  natural  aspect  ratio  for  the  fatigue  cracks  generally  fall 
in  the  range  1.3  <  a/2c  <1.6. 

The  nonsymmetric  corner  crack  results  presented  in  Figure  14  are  the 
only  data  obtained  for  that  geometry.  This  experiment  was  preliminary  in 
nature  and  was  conducted  mainly  to  determine  suitability  of  the  experimen¬ 
tal  approach  for  the  nonsymmetric  flaw  configuration.  As  indicated,  two 
corner  cracks,  intentionally  precracked  to  different  initial  sizes,  were 
located  at  opposite  sides  of  the  hole.  As  expected,  final  fracture  was 
controlled  by  the  larger  of  the  two  cracks.  In  this  particular  test,  rela¬ 
tively  few  measurements  were  obtained  during  the  transition  phase,  which 
occured  quite  rapidly  and  led  almost  immediately  to  specimen  fracture. 


2 .  Lug  Embedded  Crack  Results 

The  PMMA  lug  specimens  were  gripped  at  one  end  and  loaded  through  a 
steel  pin  placed  in  the  large  hole  in  the  other  end  as  shown  in  Figure  8. 
The  degree-of-fi t  between  the  lugs  and  the  pin  used  to  transfer  load  to 
the  lugs  was  such  that  the  pin  could  be  removed  by  hand  when  no  load  was 
applied.  Although  there  were  no  specified  tolerances  for  the  lug/pin  fit, 
all  of  the  lugs  were  machined  and  fit  by  the  same  technician. 
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The  actual  crack  lengths  for  the  embedded  flaw  lug  specimens  were 
plotted  in  the  same  manner  as  the  plate  specimens  and  are  presented  in 
Figures  16  through  Figure  24.  It  should  be  noted  that  several  of  the  lug 
tests  were  slightly  different  from  the  others  in  terms  of  the  collected 
data.  The  embedded  crack  in  specimen  PT10  (Figure  22)  grew  slightly  out 
of  plane  (approximately  15°)  so  that  the  crack  lengths  measured  in  the  2a 
direction  are  projections  of  the  actual  length.  A  corner  crack  also  exis¬ 
ted  at  the  opposite  side  of  the  hole  in  this  specimen  but  did  not  grow 
during  the  test.  The  other  exceptions  were  lug  specimens  PT14,  PT11 ,  and 
PT13  (Figures  21,  23,  and  24), where  fracture  occurred  before  both  sides  of 
the  crack  penetrated  the  lug  thickness.  Since  the  latter  specimens  con¬ 
tained  the  smallest  wall  thicknesses  (W/D  =  1.5  and  2.0),  the  embedded 
cracks  had  a  proportionately  larger  size  upon  penetrating  the  specimen 
thickness,  and  fracture  occurred  sooner. 

Note  again  in  Figures  1 6  through  23  that  the  transition  crack  dimen¬ 
sions  cL  and  cR  initially  grow  quite  rapidly  after  penetrating  the  wall 
thickness,  slow  down  as  the  middepth  crack  dimension  c  is  approached,  and 
then  speed  up  again  prior  to  final  fracture.  This  behavior  is  similar  to 
that  seen  in  the  plate  specimens,  and  occurred  both  for  nearly  symmetric 
flaws  (when  cR  Z  c^)  and  for  significantly  skewed  crack  shapes. 

The  embedded  crack  aspect  ratio  a/c  for  the  different  lug  shapes  are 
plotted  versus  flaw  size  2a/T  in  Figures  25  to  28.  Although  different 
initial  crack  shapes  were  considered,  the  natural  aspect  ratio  upon  pene¬ 
trating  the  wall  thickness  (2a/T  =  1.0)  was  almost  always  approximately 
a/c  =  1 .5. 
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Figure  16a  Fatigue  Crack  Growth  Data  for  Luc  Specimen  PT4 
(W/D  =  3.0)  Prior  to  Transition 
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Fiours  16b  Fatigue  Crack  Growth  Data  for  Ug  Specimen  PT4 
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Fatigue  Crack  Growth  Data  for  Lug  Specimen  PT7 
(W/D  =  2.5)  Prior  to  Transition 
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Figure  19b  Fatigue  Crack  Growth  Data  for  Lug  Specimen 
PT7  (W/D  =  2.5)  During  Transition  Period 
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Fatigue  Crack  Growth  Data  for  Lug  Specimen 
PT8  ( W/D  =20)  During  Transition  Period 
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Figure  22b  Fatigue  Crack  Growth  Data  for  Lug  Specimen  PT10 
(W/D  =  1.5)  During  Transition  Period 
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Figure  25  Crack  Shape  a/c  Versus  Crack  Size  2a/T  for  W/D 
Pin-Loaded  Lug  Tests 


LUG  ASPECT  RATIOS  (W/D=2.5) 

3.500 -i 


3.000 H 


□  PT6 
O  PT7 


2.500 H 


2.000  H 


1 .500  H 


i .ooo  H 


n 

m 


% 


□ 

□ 

m 


□ 


□ 


□ 

□ 


o 

_  oo 

© 

(3D 


□ 


s 


.500  H 


.ooo  H — 
.000 


T 


.200 

2fl/T 


.400 


1 

.600 


1 - 1 

.800  1 .000 


Figure  26  Crack  Shape  a/c  Versus  Crack  Size  2a/T  for  W/D  =  2.5 
Pin-Loaded  Lug  Tests 
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3.  Lug  Through  Crack  Test  Results 


Several  of  the  lug  through-the-thickness  crack  tests  were  continua¬ 
tions  of  the  embedded  tests  for  the  appropriate  specimens  (PT4,  PT5,  PT8). 
After  the  cracks  were  large  enough  to  be  considered  through-the-thickness 
cracks,  the  load  was  lowered,  if  necessary,  to  avoid  fracture,  and  the 
fatigue  tests  continued.  Lug  specimen  PT9  was  saw  cut  to  initiate  a  through- 
the-thickness  crack.  All  of  the  crack  lengths  were  determined  from  the 
test  photographs  using  the  five  point  average  discussed  earlier.  This  effec¬ 
tive  crack  length  was  called  c  to  be  consistant  with  the  other  tests.  The 
crack  length  versus  elapsed  load  cycle  data  are  shown  in  Figure  29  for  all 
of  the  tests.  (Note  that  the  specimen  PT5  data  are  given  in  two  plots). 

The  through-cracked  lug  results  were  used  to  compute  stress  intensity 
factors  from  the  known  fatigue  crack  growth  law  for  the  test  material. 

Crack  growth  rates  were  computed  by  a  least  squares  sliding  polynomial  method 
as  described  in  Section  2  for  the  baseline  tests.  Equation  1  was  then  used 
to  compute  the  apparent  cyclic  stress  intensity  factor  corresponding  to  the 
crack  length  for  that  particular  growth  rate.  These  data  are  given  in 
dimensionless  form  in  Figures  30  and  31.  Here  K  is  the  cyclic  stress  inten¬ 
sity  factor  compute"’  by  Equation  1  for  the  measured  crack  growth  rate,  D  is 
the  hole  diameter  in  the  lug,  Rq  is  the  outer  lug  radius  (R  =  W/2),  R..  is 
the  radius  of  the  hole  (R-  =  D/2),  c  is  the  average  through-the-thickness 
crack  length,  and  aQ  is  the  remote  cyclic  stress  defined  by 

P  -  P  . 

max  min  /0\ 

°o - TW -  (2) 

In  equation  2,  P  „  and  P  .  are  the  maximum  and  minimum  cyclic  loads,  T 
max  min 

is  the  lug  thickness  and  W  is  the  width  as  shown  in  Figure  1. 
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Figure  30  Comparison  of  Experimental  Measurements  of 

Dimensionless  Stress  Intensity  Factors  Computed 
from  Fatigue  Crack  Growth  Rates  with  Hsu  Finite 
Element  Analysis  (W/D  =  2.0) 
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31  Comparison  of  Experimental  Measurements  of 

Dimensionless  Stress  Intensity  Factors  Computed 
from  Fatigue  Crack  Growth  Rates  with  Hsu  Finite 
Element  Analysis  (W/D  =  3.0) 
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Stress  intensity  factors  reoorted  by  Hsu  [13]  for  these  two  through- 
cracked  lug  shapes  are  also  given  for  comparison  in  Figures  30  and  31. 

Note  that  these  numerical  results.,  obtained  by  the  finite  element  method, 
were  for  a  steel  pin  in  a  steel  lug.  These  solutions  generally  agree  fairly 
well  with  the  experimental  measurements,  which  involved  a  steel  pin  in  a 
PMMA  lug. 

Additional  details  of  the  method  used  to  compute  the  stress  intensity 
factors  from  the  measured  crack  lengths  are  given  in  Reference  11.  Stress 
intensity  factors  were  also  computed  from  the  measured  surface  crack 
growth  rates  in  the  plate  and  embedded  lug  tests  by  a  similar  procedure. 
Those  results  are  reported  separately  in  Reference  11. 
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SECTION  IV 


FATIGUE  CRACK  GROWTH  PREDICTIONS 

1 .  Approach 

The  goal  here  is  to  describe  an  approach  for  predicting  the 
growth  of  surface  or  corner  cracks  located  along  the  bore  of  a  hole 
in  a  wide  plate.  In  Figure  32,  the  surface  crack  is  represented  by 
a  semiellipse  whose  major  and  minor  axes  are  defined  by  the  coordinates 
(x-j,  y-j),  (x2,  y2),  and  (x^,  y2)  of  points  1,  2,  and  3.  The  surface 
dimension  2a  (measured  along  the  bore  of  the  hole)  and  depth  c  are 
given  by 


2a  =  Xg  -  x-|  (3) 

c  =  y2  (4) 

When  subjected  to  an  applied  cyclic  stress,  crack  tips  1,  2,  and 
3  will  grow.  Let  the  fatigue  crack  growth  rates  for  crack  tips  1, 

Jl  JO  JO 

2,  and  3  be  and  Note  that,  in  general,  these  crack  growth 

rates  will  have  different  values.  The  crack  growth  increments  for  a 
specified  number  of  cycles  aN  of  loading  are,  however,  related. 

Assume,  for  example,  that  the  depth  c  extends  a  small  amount  a2  during 
the  interval  AN.  Then,  by  definition,  ^  ,  which  gives 

AN  =  |§  (5) 

dN 

Now  the  growth  a!  of  crack  tip  1  is  given  by 

Al  =  ^  •  AN  (6) 
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Coordinate  Systei 
Surface  Crack  at 
Analysis  Scheme 


I 


Likewise,  the  extensions  of  tips  2  and  3  are  given  by 


a2  =  assumed  small  increment 


A3  =  §  •  AN 


(7) 

(8) 


Now,  if  the  cyclic  stress  intensity  factors  aK^ ,  aK^ 


aK, 


are  known  at  crack  tips  1,  2,  and  3,  the  appropriate  fatigue  crack 
growth  model  (Equation  1  for  the  present  case)  can  be  used  to  compute 
the  crack  growth  rates  and  which  in  turn  give  the  crack 

growth  increments  a1  and  a3  from  Equations  6  and  8.  The  new  positions 
(xni,  yni)  of  the  surface  crack  tips  are  then  given  by 


xnl  =  X1 


Al 


xn2  =  y2  +  a2 


xn3  =  x3  +  A3 


(9) 


Successive  iteration  of  Equations  3  to  9  gives  the  growth  of  surface 
crack  dimensions  a  and  c  as  a  function  of  elapsed  cycles  N. 

Stress  intensity  factor  solutions  reported  by  Newman  and  Raju 
[14]  for  various  surface  and  corner  crack  geometries  were  used  here 
to  compute  aK^ ,  Al^,  and  AK^.  The  Newman-Raju  results  consist  of 
empirical  equations  fit  to  finite  element  K  solutions,  and  are 
readily  programmed  for  computer  use. 


if.  .  * '"k' ».  1 • 


1 
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2. 


Results  of  Predictions 


Crack  growth  predictions  were  obtained  for  the  plate  specimens 
by  the  two  degree  of  freedom  model  described  in  the  previous  subsection 
and  are  presented  here.  Recall  that  individual  crack  tips  are 
allowed  to  grow  at  independent  rates,  so  that  the  crack  shapes 
develop  "naturally”. 

Predictions  for  plate  tests  T1 ,  T2,  and  T3  are  compared  with 
experimental  results  in  Figures  33  to  35.  Note  that  these  graphs 
are  duplicates  of  Figures  11  to  13,  except  that  the  crack  growth 
predictions  have  been  added  as  the  solid  lines.  Prior  to  penetration 
of  a  free  face  by  the  surface  crack,  calculations  are  given  for  both 
the  a  and  c  dimensions  of  the  surface  flaw  (Figures  33a,  34a,  and 
35a).  Note  that,  in  general,  the  crack  growth  is  predicted  quite 
well.  (To  avoid  extrapolating  Equation  1  beyond  its  limits,  the  ini¬ 
tial  portions  of  some  experiments  were  omitted  for  predictive  purposes). 

When  the  model  predicted  free  surface  penetration,  it  was 
assumed  the  surface  crack  instantly  changed  to  a  uniform  through- 
the-thickness  flaw  of  length  c.  The  Bowie  [15]  stress  intensity 
factor  solution  was  then  used  to  analyze  the  through-cracked  hole, 
and  the  crack  growth  predictions  continued.  The  results  for  the 
assumed  through-crack  are  compared  with  the  actual  transitioning 
flaw  shapes  in  Figures  33b,  34b,  and  35b.  Note  that  although 
the  assumed  through-crack  grows  faster  than  the  actual  flaw,  the 
total  specimen  life  is  predicted  quite  well. 
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T2 


Figure  34a  Comparison  of  Predicted  and  Actual  Crack  Growth 
Prior  to  Penetration  for  Plate  Specimen  T2 
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Figure  35a  Comparison  of  Predicted  and  Actual  Crack  Growth  Prior 
to  Penetration  for  Plate  Specimen  T3 


LENGTH 


Snow  [5]  has  presented  experimental  results  for  corner  cracked 
hole  specimens  made  from  the  same  sheet  of  PMMA  studied  here.  The 
predictive  model  was  modified  to  consider  the  corner  cracked  geometry 
and  used  to  predict  lives  for  Snow's  tests.  Those  results  are  given 
in  Figures  36  to  40.  As  before,  the  open  symbols  represent  the 
measured  crack  dimensions,  while  the  solid  lines  are  the  predicted 
results. 

Note  that  the  actual  crack  shapes  were  "pinched  in" at  the 
free  surface  so  that  crack  dimension  c  was  somewhat  less  than  length 
c*,  measured  slightly  below  the  surface  (see  sketch  on  Figure  36). 

For  analysis  purposes,  the  corner  crack  was  assumed  to  be  a  quarter 
ellipse  with  dimensions  a  and  c,  and  the  additional  crack  growth 
in  the  specimen  interior  was  ignored. 

Note  that  the  corner  crack  growth  curves  are  predicted  fairly 
well,  although  agreement  is  not  quite  as  close  as  obtained  previously 
for  the  embedded  crack  plate  results.  Nevertheless,  the  multi -degree 
of  freedom  crack  growth  model,  based  on  the  Newman-Raju  empirical 
stress  intensity  factor  equations,  makes  an  adequate  prediction 
for  crack  length  and  shape  as  a  function  of  applied  cycles. 

An  attempt  was  made  to  apply  the  crack  growth  prediction 
scheme  to  the  lug  results  described  in  the  last  section.  The  difficulty 
lies  in  the  fact  that  stress  intensity  factors  are  not  available 
in  a  general  form  for  pin-loaded  surface  cracks.  Stress  intensity 
factors  were  estimated  by  modifying  the  Newman-Raju  results  [14] 
for  embedded  cracks  at  an  open  hole  in  a  wide  plate.  A  "correction 
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Figure  38  Comparison  of  Predicted  and  Actual  Crack  Growth  for 
Snow  Corner  Crack  Test  5 
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Figure  39  Comparison  of  Predicted  and  Actual  Crack  Growth 
for  Snow  Corner  Crack  Test  6 
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Figure  40  Comparison  of  Predicted  and  Actual  Crack  Growth 
for  Snow  Corner  Crack  Test  8 


factor"  was  obtained  by  dividing  the  stress  intensity  factor 
for  a  through-cracked  lug  [13]  by  the  Bowie  result  for  a  through- 
cracked  open  hole  [15].  It  was  hoped  that  when  the  Newman-Raju 
surface  crack  result  was  multiplied  by  the  "correction"  factor, 
one  would  obtain  a  working  estimate  for  the  surface  cracked  lug 
stress  intensity  factor.  Although  crack  growth  predictions  obtained 
with  these  approximate  stress  intensity  factors  gave  fairly  reasonable 
results  in  some  cases,  the  approach  was  not  judged  successful  enough 
for  general  use. 
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SECTION  V 


DIGITIZED  CRACK  SHAPE  MEASUREMENTS 

The  objective  of  this  section  is  to  present  detailed  measurements 
of  crack  shape  changes  during  the  fatigue  testing.  Digitized  measure¬ 
ments  of  the  fatigue  crack  profiles  are  given  for  Snow's  corner  crack 
tests  [5]  and  for  one  of  the  current  embedded  surface  cracked  plates 
(specimen  T2).  It  is  believed  that  this  detailed  crack  shape  infor¬ 
mation  could  be  of  interest  to  those  who  model  surface  and  corner  cracks 
by  analytical  or  numerical  procedures. 

Since  the  crack  shapes  were  recorded  on  film  as  a  function  of  elapsed 
cycles,  it  was  possible  to  make  detailed  measurements  of  the  crack  pro¬ 
files.  As  shown  schematically  in  Figures  41  and  42,  tabulated  results 
are  given  here  for  the  coordinates  of  points  along  the  perimeter  of 
corner  and  embedded  surface  cracks  located  at  the  bore  of  the  hole  in 
the  PMMA  plates.  Digitized  measurements  are  given  for  nine  points  along 
the  corner  and  transitioning  corner  cracks  (Figure  41),  while  the  co¬ 
ordinates  of  ten  embedded  surface  crack  and  transitioning  flaw  points 
are  presented  (Figure  42). 

The  corner  crack  results  were  obtained  by  re-examining  Snow's  orig¬ 
inal  filmstrips  on  a  photo  interpreter/digitizer.  The  insturment  auto¬ 
matically  punched  digitized  measurements  of  the  crack  profiles  onto 
computer  cards  for  subsequent  analysis.  The  data  reduction  system  was 
estimated  to  locate  points  along  the  crack  profile  to  an  accuray  of  +_ 
0.003  inches  (actual  size).  A  different  instrument  was  used  to  measure 
the  films  for  embedded  crack  specimen  T2  tested  here.  In  this  case 
a  Mann  digital  comparator  was  employed,  and  again  estimated  accur- 
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Schematic  Representation  of  Embedded  Surface  and  Through-Thickness 
Cracks  Showing  Location  of  Digitized  Points  Along  Crack  Front 


acy  of  +  0.003  inches  (actual  size)  was  obtained.  The  latter  instru¬ 
ment  did  not  have  the  automatic  data  recording  feature,  however,  and 
the  digitized  coordinates  had  to  be  recorded  by  hand,  and  then  entered 
into  the  computer  for  analysis.  Since  this  was  a  fairly  time  consum¬ 
ing  procedure,  only  one  of  the  embedded  surface  crack  tests  was  exam¬ 
ined  in  this  detail. 

In  all  cases,  the  photographs  were  scaled  with  respect  to  known 
distances  in  the  x  and  y  directions.  Snow's  original  measurements  for 
the  a  and  c  crack  dimensions  reported  in  Reference  5  were  based  on  a 
scaling  dimension  in  the  x  direction  only,  since  he  assumed  that  the 
viewing  mirror  was  located  at  45°  to  the  crack  plane,  and  thus  gave 
a  full  view  of  the  crack  surface.  Subsequent  work  with  the  loading 
apparatus  in  the  present  experiments  indicated  that  the  viewing  mirror 
could  sometimes  give  an  oblique  view  of  tne  crack  plane.  By  using  scaling 
dimensions  in  two  perpendicular  directions,  it  was  possible  to  deter¬ 
mine  an  angle  correction  and  to  compute  the  actual  crack  place  measure¬ 
ment.  All  crack  length  data  reported  here  have  been  corrected  in  this 
manner. 

Digitized  data  for  Snow's  [5]  corner  crack  Test  6  are  given  in 
graphical  and  tabular  form  in  Figure  43  and  Table  3  (this  information 
for  Test  6  is  reproduced  here  from  Reference  6).  Here  the  remote  stress 
varied  between  18  and  630  osi,  the  hole  diameter  D  =  0.739  in.,  the  total 
width  W  =  7.950  in.,  and  the  plate  thickness  T  =  0.698  in.  Note  from 
Figure  43  that  the  crack  size  >nd  shapes  are  well  documented  as  the  cor¬ 
ner  crack  grows  along  the  bore  of  the  hole  and  transitions  through  the 
thickness  of  the  plate.  Fatigue  crack  growth  data  obtained 
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Table  3  continued 
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Table  3  continued 
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23700  0.0  .089  .176  .253  .298  .330  .352  .369 

.670  .642  .576  .487  .400  .325  .242  .078 
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24800  0.0  .553  .609  .622  .633  .635  .629  .615 

.698  .698  .545  .462  .356  .287  .146  .060 


24900  0.0  .938  .994  1.011  1.018  1.023  1.015  1.001  .978 

.698  .698  .574  .475  .405  .293  .140  .063  0.0 


Table  3  continued 


i 


from  the  same  sheet  of  PMMA  at  the  1Hz  cyclic  frequency  employed  by  Snow  are 
given  in  Figure  44.  Similar  tabular  and  graphical  data  are  given  in  the 
Appendix  for  Snow  corner  crack  Tests  1,  3,  5,  8,  and  9.  As  indicated  previously, 
these  measurements  were  all  taken  from  Snow's  original  filmstrips  and  have  been 
corrected  for  misalignment  of  the  viewing  mirror. 

Digitized  crack  profiles  for  embedded  surface  crack  Test  T2  are  given  in 
Table  4  and  in  Figure  45.  Note  that  only  selected  crack  profiles  are  detailed 
here.  Additional  measurements  of  the  a  and  c  dimensions  (obtained  by  indepen¬ 
dent  measurement  with  the  film  strip  projector  method  described  in  Section  II) 
were  given  earlier  in  Figures  12a  and  12b.  Again  note  that  the  crack  shape 
changes  are  characterized  quite  well  by  the  digitized  profiles  shown  in  Figure 
45. 

Additional  discussion  of  the  behavior  of  the  Snow  corner  crack  test  6  and 
the  present  embedded  flaw  Test  T2  are  given  in  [12].  In  particular.  Figures 
6  and  8  in  [12]  present  stress  intensity  factors  obtained  by  the  three-dimen¬ 
sional  finite  element-alternating  method  for  particular  crack  shapes.  The 
stress  intensity  factor  changes  are  correlated  with  fatigue  crack  growth  rate 
behavior  at  the  ends  of  the  crack  tips. 
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Digitized  Measurements  of  Crack  Profiles  for  Embedded  Surface  Crack  Test  T2 


SECTION  VI 


SUMMARY  AND  CONCLUSIONS 

A  series  of  fatigue  tests  have  been  conducted  with  initial  embedded 
surface  cracks  located  along  the  bore  of  a  hole  in  a  wide  plate  or  a  pin- 
loaded  lug.  The  tests  were  conducted  with  transparent  polymer  specimens 
which  allowed  in  situ  observation  of  the  crack  plane.  Crack  growth  was 
recorded  by  time  lapse  photography.  Subsequent  measurement  of  the  crack 
photographs  gave  crack  size  and  shape  changes  as  a  function  of  elapsed 
cycles. 

A  multidegree  of  freedom  fracture  mechanics  model  was  developed  to 
predict  the  growth  of  embedded  surface  and  corner  cracks  located  at  open 
holes  in  plates  loaded  in  remote  tension.  The  analysis  employed  empirical 
stress  intensity  factor  solutions  reported  by  Newman  and  Raju  [14]  for 
surface  and  corner  cracks  located  at  open  holes.  Crack  shape  was  a  free 
parameter  in  the  predictive  model,  and  flaw  shapes  were  allowed  to  develop 
naturally.  Upon  penetration  of  a  free  surface,  the  surface  and  corner 
cracks  were  assumed  to  instantly  transition  into  uniform  through-the- 
thickness  flaws. 

Predictions  by  the  multidegree  of  freedom  model  gave  excellent  agree¬ 
ment  with  the  actual  growth  of  embedded  surface  cracks  at  open  holes.  The 
initial  growth  of  the  width  2a  and  depth  c  of  the  embedded  surface  cracks 
agreed  very  well  with  the  predictions.  The  assumed  transition  into  an 
instant  through-crack  gave  a  conservative,  but  good  estimate  of  total 
specimen  life.  The  predictive  model  was  also  used  to  analyze  corner 
cracked  hole  experiments  reported  earlier  by  Snow  [5].  Although  the 
corner  crack  predictions  generally  gave  shorter  lives  than  observed  experi¬ 
mentally,  and  did  not  show  as  close  agreement  as  the  embedded  surface 
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cracks,  the  corner  crack  analysis  still  gave  reasonable  predictions  for 
total  specimen  life.  Thus,  it  is  concluded  that  the  Newman  and  Raju  [14] 
stress  intensity  factor  solutions  for  surface  and  corner  cracked  holes  are 
quite  useful  for  predicting  the  shape  and  size  of  surface  and  corner  cracks 
at  holes.  An  attempt  to  modify  the  Newman-Raju  open  hole  solutions  with  a 
"correction  factor"  in  order  to  apply  them  to  the  pin-loaded  tests  met  with 
only  limited  success.  The  "correction  factor"  approach  is  not  recommended 
for  general  use. 

The  Snow  corner  crack  films  and  one  of  the  present  set  of  embedded 
surface  crack  photos  were  measured  in  a  manner  which  gave  a  detailed  record 
of  the  change  in  crack  shape  with  fatigue  life.  The  crack  shapes  were 
digitized,  and  coordinates  of  points  along  the  crack  border  are  presented 
in  both  graphical  and  tabular  format.  These  data  provide  a  detailed  record 
of  naturally  occurring  fatigue  crack  shapes  and  may  be  of  interest  to  those 
who  model  surface  or  corner  cracks  by  numerical  techniques. 

The  transition  of  corner  and  surface  cracks  into  uniform  through-the- 
thickness  flaws  was  examined  in  detail.  Both  three-dimensional  stress 
intensity  factor  solutions  and  the  fatigue  crack  growth  results  indicate 
that  initial  corner  or  embedded  surface  cracks  try  to  grow  into  a  stable 
through-the-thickness  configuration.  The  stress  intensity  factor  varies 
along  the  crack  perimeter  so  that  the  maximum  K  occurs  at  the  trailing, 
rather  than  the  leading  point  of  crack  advance.  Moreover,  at  the  trailing 
point,  where  the  crack  penetrates  a  free  surface,  K  varies  significantly 
with  crack  advance. 

The  stress  intensity  factor  reaches  a  large  value  initially  after 
free  surface  penetration,  but  then  decreases  locally  at  that  point  as  the 
fre'  face  crack  length  grows  to  a  uniform  through-the-thickness  shape. 
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The  stress  intensity  factor  at  the  leading  point  of  maximum  crack  advance 
usually  has  a  smaller  magnitude  than  at  the  trailing  point,  and  is 
apparently  unaffected  by  the  large  changes  occurring  at  the  trailing 
position  during  transition.  The  leading  edge  stress  intensity  factor  may 
often  be  approximated  during  the  transition  period  by  the  two-dimensional 
analysis  for  a  uniform  through-the-thickness  crack  whose  length  equals  the 
distance  of  maximum  crack  advance. 

The  measured  fatigue  crack  growth  rates  agree  with  the  computed  stress 
intensity  factors,  as  the  crack  perimeter  advances  locally  at  rates  corres¬ 
ponding  to  the  K  variation  along  the  flaw  perimeter.  In  particular,  the 
free  face  crack  dimension  iirmediately  grows  quite  rapidly  after  penetration, 
but  then  decreases  to  the  rate  seen  by  the  point  of  maximum  crack  advance. 


APPENDIX 


DIGITIZED  MEASUREMENTS  OF 
SNOW  CORNER  CRACK  TESTS 

This  appendix  presents  digitized  measurements  of  corner  cracked  hole 
experiments  conducted  by  Snow  [5].  The  original  filmstrips  were  remeasured, 
and  the  results  given  here  in  graphical  and  tabular  form.  The  format  is  the 
same  as  that  described  in  Section  V  for  Test  6.  Results  are  given  for 
Snow's  Test  Numbers  1,  3,  5,  8,  and  9.  Specimen  dimensions  are  given  in 
reference  [12].  In  all  cases,  the  remotely  applied  tension  force  varied 
between  100  and  3500  lbs.  so  that  the  applied  cyclic  force  was  3400  lbs. 

At  the  beginning  of  each  table,  NUM  refers  to  the  total  number  of  crack 
shapes  digitized,  2R  is  the  hole  diameter,  W  is  the  plate  width,  T  is  the 
plate  thickness,  and  LOAD  is  the  applied  cyclic  force. 
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Table  5  Digitized  Data  for  Snow  Test  1 


NUM  =  S5  2R  =  .747  W  =  7.380 
T  =  .729  LOAD  =  3400. 
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0 
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Table  5  continued 
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.150 

.159 

.153 

.140 

.177 
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0 
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0 
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Table  5  continued 


CYCLES  = 

I  = 

1S500. 

1 
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9 
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.070 

.109 

.147 

.175 

.193 

.201 

.193 

.181 

Y(I)  = 

.344 

.321 

.253 

.233 
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.141 

.089 

.030 

0 
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I  = 
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1 

2 

3 

4 

5 
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3 
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.053 

.105 

.138 
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.199 
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1 
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3 

9 
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0 
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1 

2 

3 

4 
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3 
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0 

.059 

.103 

.144 

.180 

.202 

.211 
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.332 

.348 

.312 

.239 

.211 

.151 
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.033 
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CYCLES  = 
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3 
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5 
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7 

3 
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0 
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.209 

.194 
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.303 
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.223 

.178 

.102 

.033 

0 

CYCLES  = 

I  = 

13500. 

1 

2 

3 

4 

e 
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7 

3 
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.073 
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.133 
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.220 

.212 
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.303 

.252 

.202 
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.085 

.032 

0 

CYCLES  = 

I  = 

13700. 

1 

2 

3 

4 

5 

3 

7 

8 

9 
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.110 
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.349 

.304 

.257 

.222 

.174 

.113 

.048 

0 

CYCLES  = 

I  = 

17100. 

1  2 

3 

4 

5 

3 

7 

8 

9 

Xd)  = 

0 

.075 

.125 

.  135 

.194 

.222 

.229 

.220 

.208 

Yd)  = 

.387 

.337 

•  WlS 

.271 

.224 

.130 

.083 

.028 

0 

CYCLES  =  17300. 

I  =  1  2 

4 

c 

S 

7 

3 

9 

Xd)  =  0  .074 

•  127 

.  135 

.207 

.227 

.234 

.225 

.213 

Yd)  =  .388  .333 

•  Cc3 

.273 

.203 

.  145 

.079 

.027 

0 
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Table  5  continued 


CYCLES 


17500. 


I  = 

1 

2 

3 

4 

5 

6 

7 

3 

9 

X(I)  = 

0 

.084 

.117 

.184 

.199 

.227 

.237 

.234 

.217 

yen  = 

.339 

.387 

.347 

.292 

.236 

.182 

.110 

.046 

0 

CYCLES  = 

17700. 

I  = 

1 

2 

3 

4 

5 

6 

7 

8 

9 

xcn  = 

0 

.085 

.133 

.172 

.205 

.228 

.244 

.240 

.226 

Y(I)  * 

.420 

.392 

.354 

.302 

.240 

.184 

.105 

.044 

0 

CYCLES  = 

17S00 . 

I  = 

1 

2 

3 

4 

5 

6 

7 

3 

9 

X  C I )  = 

0 

.092 

.140 

.179 

.208 

.234 

.248 

.242 

.224 

YCI)  = 

.428 

.  3S9 

.358 

.298 

.244 

.178 

.100 

.046 

0 

CYCLES  = 

18100. 

I  = 

1 

2 

3 

4 

5 

6 

7 

8 

9 

xcn  = 

0 

.088 

.150 

.183 

.224 

.245 

.251 

.245 

.229 

yen  = 

.433 

.408 

.349 

.2S8 

.218 

.155 

.098 

.037 

0 

CYCLES  = 

18300. 

I  = 

1 

2 

3 

4 

5 

6 

7 

8 

9 

xcn  = 

0 

.088 

.151 

.190 

.217 

.234 

.253 

.251 

.235 

Y(I)  a 

.440 

.412 

.  358 

.2S4 

.243 

.202 

.118 

.052 

0 

CYCLES  = 

18500. 

T  s 

1 

2 

3 

4 

5 

6 

7 

3 

9 

xcn  = 

0 

.090 

.174 

.207 

.235 

.243 

.259 

.258 

.238 

YCI)  = 

.447 

.421 

.335 

.274 

.215 

.177 

.118 

.052 

0 

CYCLES  = 

18700. 

I  a 

1 

2 

3 

4 

5 

6 

7 

8 

9 

xcn  a 

0 

.092 

.144 

.ISO 

.224 

.248 

.287 

.262 

.246 

YCI)  a 

.458 

.430 

.378 

.321 

.250 

.195 

.108 

.050 

0 

CYCLES  = 

1SS00. 

I  = 

1 

2 

3 

4 

5 

6 

7 

8 

9 

xcn  = 

0 

.083 

.144 

.197 

.230 

.250 

.289 

.  253 

.249 

YCI)  a 

.453 

.442 

.387 

.315 

.  248 

.197 

.111 

.052 

0 

CYCLES  = 

13100. 

I  -a 

1 

2 

3 

4 

5 

6 

7 

3 

9 

xcn  = 

0 

.078 

.131 

.184 

.227 

.253 

.273 

.287 

.254 

YCI)  a 

.470 

.448 

.408 

.  348 

.  284 

.201 

.100 

.042 

0 

CYCLES  = 

1SEC0 . 

T  = 

2 

4 

c 

_S 

7 

3 

9 

xcb  a 

c 

.CSC 

.138 

.183 

.223 

.  254 

.279 

.272 

.259 

YCI)  = 

.472 

.410 

•  LC  < 

.  273 

.215 

.  102 

.  033 

0 
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Table  5  continued 


CYCLES 
I  = 

=  19500. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

X(  1 5  = 

0 

•  063 

.  142 

.212 

.243 

.276 

.286 

.276 

.267 

Y(I)  = 

.463 

.474 

.421 

.323 

.246 

.175 

.096 

.030 

0 

CYCLES 
I  = 

=  19700. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

XCI)  = 

0 

.OSS 

.169 

.243 

.267 

.283 

.290 

.282 

.270 

YC 1 3  = 

.488 

.466 

.  3S4 

.268 

.212 

.151 

.095 

.043 

0 

CYCLES 
I  = 

=  20100. 
1 

2 

3 

4 

5 

6 

7 

3 

9 

XCI)  = 

0 

.113 

.ISO 

.243 

.260 

.235 

.301 

.290 

.280 

YC  1 3  = 

.510 

.472 

.368 

.253 

.210 

.168 

.092 

.025 

0 

CYCLES 
I  = 

=  20300. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

XCI)  = 

0 

.122 

.218 

.257 

.282 

.302 

.307 

.303 

.288 

YC  I J  = 

.520 

.477 

■  366 

.  282 

.223 

.156 

.095 

.039 

0 

CYCLES  =  20500. 

I  =  1 

2 

3 

4 

5 

6 

7 

8 

9 

XCI)  =  0 

.109 

.165 

.231 

.278 

.297 

.316 

.307 

.296 

YC I 3  =  .534 

.453 

.425 

.356 

.252 

.201 

.100 

.027 

0 

CYCLES  =  20700. 

I  =  1 

2 

r> 

4 

5 

6 

7 

8 

9 

XCI)  =  0 

.121 

.204 

.255 

.251 

.313 

.322 

.314 

.300 

-C 

n 

II 

l/l 

CJ1 

.505 

.416 

.318 

.236 

.ISO 

.106 

.035 

0 

CYCLES  =  20SOO. 

I  =  1 

2 

2 

4 

5 

6 

7 

3 

9 

XCI)  =  0 

.122 

.205 

.2SS 

.237 

.320 

.331 

.32S 

.310 

YC I )  =  .553 

.515 

.430 

.315 

.246 

.174 

.  0S6 

.043 

0 

CYCLES  =  21100. 

I  =  1 

2 

3 

4 

5 

6 

7 

3 

9 

XCI)  =  0 

.129 

.204 

.  286 

.237 

.321 

.336 

.323 

.312 

YC I )  =  .565 

.543 

.462 

.353 

.281 

.197 

.106 

.032 

.002 

CYCLES 
I  = 

=  21300. 
1 

2 

3 

4 

5 

6 

7 

3 

9 

XCI)  = 

0 

.134 

.215 

.231 

.306 

.330 

.342 

.332 

.323 

YC  I )  = 

.552 

.544 

.454 

.342 

.267 

.186 

.083 

.023 

0 

CYCLES 

--  21500. 

1 

2 

4 

c 

6 

7 

3 

9 

XCI)  = 

C 

.125 

•  21S 

.253 

.325 

.341 

.  350 

.341 

.328 

YC  I)  = 

.603 

.533 

.  4£S 

.  331 

.240 

.174 

.106 

.026 

0 
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Table  5  continued 


CYCLES  = 

I  = 

21700. 

1 

2 

3 

4 

5 

6 

7 

3 

9 

Xd)  a 

0 

.126 

.247 

.311 

.337 

.351 

.363 

.351 

.342 

Yd)  = 

.614 

.573 

.454 

.322 

.242 

.187 

.105 

.025 

0 

CYCLES  = 

I  = 

21300. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Xd)  3 

0 

.138 

.242 

.306 

.336 

.358 

.371 

.365 

.349 

Yen  = 

.632 

.583 

.483 

.357 

.280 

.193 

.118 

.040 

0 

CYCLES  * 

I  = 

22100. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

X  C 1 5  = 

0 

.139 

.236 

.318 

.352 

.376 

.384 

.375 

.366 

Y(I)  = 

.652 

.609 

.514 

.363 

.274 

.183 

.091 

.035 

0 

CYCLES  = 
1  = 

22300. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

XCI 5  = 

0 

.165 

.245 

.314 

.354 

.380 

.333 

.383 

.369 

Yd)  = 

.670 

.608 

.525 

.402 

.290 

.198 

.111 

.031 

0 

CYCLES  = 

I  = 

22500. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Xd)  = 

0 

.151 

.256 

.313 

.362 

.386 

.405 

.396 

.383 

Yd)  * 

.683 

.636 

.540 

.441 

.304 

.224 

.114 

.031 

0 

CYCLES  = 

I  = 

22700. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

X(I)  = 

0 

.146 

.247 

.320 

.374 

.403 

.425 

.414 

.397 

Yd)  = 

.725 

.680 

.605 

.498 

.  366 

.257 

.127 

.040 

0 

CYCLES  = 

I  = 

22S00. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

XCI)  = 

.352 

.404 

.438 

.449 

.457 

.461 

.466 

.447 

.412 

Yd)  = 

.728 

.657 

.548 

.488 

.410 

.306 

.197 

.069 

0 
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7  Digitized  Crack  Profiles  for  Snow  Corner 
Crack  Test  3 


Table  6  Digitized  Data  for  Snow  Test  3 


I  = 
XU)  = 
yen  = 

1 

0 

.162  . 

2 

033 

162 

3 

.064 

.157 

4 

.097 

.144 

3 

.117 

.132 

.140 

.107 

.161 

.075 

.167 

.OSO 

.169 

0 

CYCLES  = 
i  = 
XU)  * 
Y(I )  = 

2000. 

1 

0 

•  186  . 

2 

034 

ISO 

3 

.071 

.184 

4 

.099 

.172 

S 

.124 

.151 

6 

.151 

.114 

7 

.170 

.082 

8 

.170 

.046 

9 

.170 

0 

CYCLES  = 
I  = 

XC I ?  = 
YCI)  * 

3000. 

1 

0 

.199 

2 

.050 

.207 

3 

.033 
.  1S5 

4 

.122 

.170 

5 

.144 

.150 

6 

.165 

.118 

7 

.178 

.087 

3 

.182 

.049 

9 

.180 

0 

CYCLES  = 
I  => 
XU)  = 
YU  )  = 

4000. 

1 

0 

.221 

2 

.  036 
.224 

3 

.087 

.213 

4 

.126 

.181 

5 

.162 

.140 

6 

.182 

.104 

7 

.188 

.060 

3 

.187 

.026 

9 

.181 

0 

CYCLES  = 
I  = 
X(I)  = 
YU)  * 

=  5C00. 

1 

0 

.245 

2 

.048 

.247 

3 

•  0S0 
.236 

4 

.125 

.211 

S 

.162 

.163 

6 

.185 

.129 

7 

.183 

.071 

3 

.197 

.026 

9 

.191 

0 

CYCLES  : 

I  = 
X(I)  = 
YU)  = 

=  sooo. 
1 

0 

.261 

2 

.043 

.263 

.OSS 

.247 

4 

.135 

.221 

5 

.162 

.1SS 

6 

.183 

.140 

7 

.207 

.084 

3 

.206 

.025 

9 

.200 

0 

CYCLES 

I  = 

xci)  = 

Y(I)  = 

=  7000 

1 

0 

.250 

• 

2 

.054 

.283 

3 

.057 

.275 

4 

.144 

.242 

5 

.164 

.166 

6 

.205 

.139 

7 

.218 

.091 

8 

.219 

.046 

9 

.208 

0 

CYCLES 

XU)  = 
YU)  = 

=  8000 
1 

.000 

.318 

2 

.050 

.318 

3 

.112 

.297 

4 

.  163 
.255 

5 

.187 

.220 

6 

.212 

.175 

7 

.228 

.113 

3 

.233 

.063 

9 

.219 

0 

CYCLES 

I  = 

XU)  = 
YU)  = 

=  9000 

1 

0 

.347 

2 

.046 

.350 

3 

.103 

.332 

4 

.ISO 

.2SS 

5 

.ISO 

.250 

6 

.216 

.202 

7 

.241 

.140 

8 

.250 

.085 

9 

.235 

0 

CYCLES 

XCI)  = 
YU)  = 

=  SEC  C 

l 

2 

.032 

.362 

"3 

.107 

.335 

4 

.173 

.234 

5 

.206 

.240 

_S 

.251 

.136 

7 

.  253 
.073 

3 

.233 
.  032 

9 

.250 

0 

m 


Table  6  continued 


CYCLES  = 

I  = 

10000. 

1  2 

3 

4 

5 

6 

7 

8 

9 

Xd)  = 

0 

.OBI 

.118 

.1B4 

.217 

.243 

.268 

.271 

.259 

Yd)  = 

.37S 

.373 

.350 

.314 

.248 

.134 

.119 

.072 

0 

CYCLES  = 

I  = 

10500. 

1 

2 

3 

4 

S 

6 

7 

3 

9 

X(I)  = 

0 

.045 

.110 

.158 

.201 

.243 

.264 

.278 

.267 

YCI)  = 

.335 

.  3S8 

.376 

.340 

.255 

.219 

.162 

.087 

0 

CYCLES  = 

I  = 

11000. 

1 

2 

3 

4 

5 

6 

7 

3 

9 

X(I)  = 

0 

.051 

.123 

.204 

.253 

.276 

.292 

.239 

.277 

Yen  = 

.422 

.425 

.  3S8 

.318 

.240 

.169 

!  103 

.040 

0 

CYCLES  = 

I  = 

11500. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

X(I)  = 

0 

.043 

.157 

.219 

.254 

.287 

.299 

.300 

.290 

Y(I)  = 

.441 

.443 

.389 

.315 

.259 

.165 

.105 

.070 

0 

CYCLES  * 

I  = 

11750. 

i 

2 

3 

4 

5 

5 

7 

8 

9 

xcn  = 

0 

.047 

.148 

.158 

.242 

.272 

.294 

.311 

.298 

Yd)  = 

.452 

.456 

.410 

.  383 

.302 

.245 

.180 

.073 

0 

CYCLES  = 

I  = 

12000. 

1 

2 

O 

4 

5 

5 

7 

3 

9 

XCI)  = 

.003 

.063 

.185 

.243 

.275 

.301 

.314 

.317 

.307 

YCI)  = 

.474 

.472 

.3SS 

.322 

.  251 

.187 

.123 

.079 

0 

CYCLES  = 

T  = 

12250. 

1 

2 

3 

4 

— 5 

6 

7 

8 

9 

Xd)  = 

.000 

.086 

•  ISO 

.236 

.232 

.298 

.  >20 

.324 

.316 

Yd)  = 

.434 

.484 

.413 

.359 

.306 

.228 

.143 

.080 

0 

CYCLES  = 

I  = 

12500. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

X(I)  = 

0 

.090 

.195 

.237 

.282 

.312 

.330 

.334 

.327 

Yd)  = 

.497 

.489 

.407 

.361 

.280 

.212 

.136 

.069 

0 

CYCLES  * 

T  s 

12750. 

1 

2 

3 

4 

5 

6 

7 

3 

9 

xd)  = 

0 

.  OSS 

.132 

.  cSS 

.284 

.318 

.338 

.342 

.333 

YCI)  = 

.506 

.510 

.477 

.382 

.283 

.205 

.  129 

.067 

0 

CYCLES  * 

13CC0. 

1 

2 

n 

w 

<. 

c 

S 

7 

3 

9 

xd)  = 

0 

.083 

.  *4S 

•  src 

•  283 

.  322 

.343 

.333 

.345 

Yd)  = 

.523 

.525 

.42  c 

.  ZZz 

.  2i? 

.  223 

.144 

.075 

0 

CYCLES  = 

I  = 

13200. 

1 

2 

3 

4 

5 

6 

7 

3 

9 

XCI)  = 

.001 

.090 

•  162 

.248 

.288 

.324 

.352 

.363 

.351 

Yd)  = 

.543 

.532 

.475 

.387 

.  326 

.239 

.150 

.062 

0 
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Table  7  Digitized  Data  for  Snow  Test  5 


MUM  a  25  2f?  =  .742  W  *  7.S50 

T  =  . 720  LOAD  *  3400. 


CYCLES  =  2000. 

I  =  1  2  3 
XCI)  =  0  . 022  .052 
YCI)  =  . 1S8  . 1SS  .1E2 


4 

.084 

.144 


S 

.103 
.  126 


B 

.119 

.109 


7 

.127 

.090 


3 

.132 

.OSS 


9 

.132 

0 


CYCLES  =  3000. 

I  =  1  2 
XC I J  =  0  .044 
Y ( I )  =  .181  .178 


3  4  5 
.071  .095  .112 
. 1G5  . 145  . 127 


S 

.125 

.107 


7 

.136 

.077 


8 

.142 

.037 


9 

.142 

0 


CYCLES  =  40CC. 

I  =  1  2  3  4 
X(T)  =  0  .040  .071  .ICO 
Y(b  =  -188  .185  .172  .150 


S 

.116 
.  129 


6 

.133 

.101 


7 

.141 

.075 


8 

.147 

.035 


9 

.146 

0 


CYCLES  =  5000 

I  =  1 

XC 1 5  =  0 

Yen  =  .187 


2  3  4 

032  .076  .109 

187  . 169  .138 


5  6  7 

128  . 145  . 143 

109  .075  .049 


3  9 

.148  .150 

.027  0 


CYCLES  = 

X(I)  * 
Yd)  = 

6000. 

1  2 

0  .041 

.207  .204 

.074 

.193 

4 

.104 

.172 

CYCLES  = 

T  = 

xcn  = 

YC 1 5  = 

7000. 

1  2 

0  .  Ow3 

.216  .212 

3 

.070 

.203 

4 

.OSS 

.187 

CYCLES  = 

I  = 
«!)  = 
Yen  = 

8000. 

1  2 

0  .043 

.22S  .223 

•J 

.  076 

.216 

4 

.104 

.199 

CYCLES  = 

T  = 

xcn  = 
y  i)  = 

9000. 

1  2 

0  .047 

.239  .241 

.077 

.231 

4 

.114 
.  2C6 

5 

6 

7 

3 

9 

127 

.142 

.153 

.155 

.155 

i  37 

.106 

.069 

.040 

0 

5 

125 

154 

6 

.141 

.124 

7 

.153 

.091 

3 

.ISO 
•  030 

9 

.160 

0 

5 

127 

170 

6 

.147 

.140 

7 

.161 

.099 

8 

.163 

.066 

9 

.167 

0 

5 

S 

7 

3 

9 

140 

.159 

.174 

.179 

.175 

171 

.139 

.083 

.045 

0 

CYCLES  = 

T  a 

10000. 

1  2 

n 

w 

xcb  = 

.000 

.049 

.CS7 

YC)  = 

.252 

.254 

.  e-5 

4 

.127 

.208 


6 

.163 
.  132 


7 

.183 

.080 


3 

.133 

.043 


9 

.181 

0 


CYCLES  =  11000. 

I  =  1  2 
XC I )  =  0  .040 
Y(I)  *  .258  .253 


2  4  5 
. CS7  .136  .ISO 
.244  .214  .178 


S  7  J 
.176  .189  .193 
.142  .103  .055 


114 


Table  7  continued 


CYCLES  = 
I  = 

12000. 

1 

2 

3 

4 

5 

S 

7 

8 

8 

XCI 3  = 

0 

.043 

.103 

.140 

.168 

.187 

.205 

.210 

.205 

Y(I)  = 

.279 

.283 

.284 

.232 

.135 

.  156 

.098 

.054 

0 

CYCLES  = 

I  = 

13000. 

1 

2 

3 

4 

5 

6 

7 

3 

9 

XC 1 3  = 

0 

.041 

.086 

.133 

.176 

.202 

.216 

.221 

.216 

YCI)  = 

.305 

.308 

.2S5 

.263 

.217 

.162 

.112 

.067 

0 

CYCLES  = 

I  = 

14000. 

1 

2 

3 

4 

5 

6 

7 

8 

3 

XCI)  = 

0 

.055 

.052 

.134 

.168 

.138 

.222 

.233 

.226 

Yd)  = 

.336 

.337 

.322 

.232 

.235 

.208 

.145 

.067 

0 

CYCLES  = 

I  = 

15000. 

1 

2 

3 

4 

5 

6 

7 

3 

9 

XCI)  = 

0 

.052 

.106 

.149 

.168 

.217 

.238 

.255 

.249 

Yd)  = 

.368 

.368 

.350 

.314 

.272 

.222 

.168 

.078 

0 

CYCLES  = 

I  = 

16000. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

XCI)  = 

0 

.052 

.104 

.186 

.228 

.247 

.262 

.273 

.267 

Y(I)  = 

.404 

.406 

.382 

.313 

.244 

.206 

.155 

.074 

0 

CYCLES  = 

I  = 

17000. 

1 

2 

3 

4 

5 

6 

7 

3 

3 

XCI)  = 

0 

.066 

.128 

.183 

.229 

.262 

.283 

.289 

.280 

Yd )  = 

.458 

.454 

.425 

.367 

.305 

.227 

.145 

.063 

0 

CYCLES  = 

I  = 

180C0. 

1  2 

3 

4 

5 

6 

7 

8 

9 

XCI)  = 

0 

.050 

.116 

.173 

.220 

.275 

.303 

.317 

.308 

YCI)  = 

.537 

.531 

.514 

.465 

.403 

.304 

.204 

.039 

0 

CYCLES  =  18250. 

I  =  1  2  3  4 
X( 1 5  =  C  .OSS  .153  .201 
YCI )  =  .547  .534  .487  .433 


CYCLES  =  18500. 

I  =  1  2  3  4 
XCI)  =  0  .058  .113  .175 
YCI)  =  .573  . 5E4  .541  .487 


CYCLES  =  187S0. 

I  =  1  2  3  4 
XC 13  =  .002  .CSC  .132  .210 
Yd)  =  .ScS  .574  . 52S  .488 


5 

6 

7 

3 

9 

.25S 

.285 

.310 

.324 

.314 

• 

.290 

.197 

.090 

0 

5 

6 

7 

3 

9 

.228 

.268 

.303 

.330 

.320 

.425 

.  350 

.254 

.103 

0 

5 

3 

r 

3 

9 

•  332 

•  -  4 

.332 

.342 

.330 

.  2S4 

.172 

.095 

0 

115 


Table  7  continued 


CYCLES  = 

I  = 

1S000. 

1 

2 

3 

4 

xcn  = 

0 

.073 

.140 

.152 

Yd)  = 

.SCO 

.533 

.535 

.505 

CYCLES  = 

I  = 

15250. 

1  2 

3 

4 

X(I)  = 

0 

.107 

.  17G 

.230 

Yd)  = 

.549 

.629 

.570 

.506 

CYCLES  = 

I  = 

19500. 

1 

2 

3 

4 

Xd)  = 

0 

.093 

.171 

.234 

YCI)  = 

•  632 

.650 

.534 

.518 

CYCLES  = 

I  = 

13750. 

1 

2 

3 

4 

xd)  = 

0 

.032 

.142 

.223 

Y(I)  = 

.698 

.  686 

.  653 

.572 

CYCLES  = 

I  = 

20000. 

1 

2 

3 

4 

Xd)  = 

0 

.139 

.263 

.306 

Yd)  = 

.715 

.715 

.616 

.517 

S 

6 

7 

3 

9 

257 

.305 

.330 

.350 

.339 

414 

.311 

.230 

.097 

0 

5 

6 

7 

8 

9 

28 1 

.314 

.341 

.356 

.345 

421 

.344 

.257 

.120 

0 

5 

6 

7 

8 

9 

260 

.332 

.358 

.367 

.357 

442 

.316 

.211 

.116 

0 

5 

6 

7 

8 

9 

277 

.329 

.362 

.332 

.373 

437 

.374 

.270 

.120 

0 

5 

6 

7 

3 

9 

330 

.362 

.337 

.398 

.386 

450 

.348 

.223 

.112 

0 

116 


C  YCLES 
26 OOO- 


-25500- 
252  5  0- 
2looo- 
2  4  7SO- 
24  Coo- 
24  300- 

2*ooO- 

2JSO0- 

2?ooO- 

Z1S0O- 


22000 


Table  8  Digitized  Data  for  Snow  Test  8 

NUM  =31  2R  =  .743  U  =  7.S80 
T  =  .704  LOAD  =  3400. 


CYCLES  = 

I  = 

XU)  = 

Y(I)  = 

3000. 

1 

0 

167  . 

2 

036 

164 

3 

.057 

.146 

4 

.082 

.120 

5 

.104 

.083 

6 

.112 

.062 

7 

.124 

.043 

8 

.127 

.014 

9 

.127 

0 

CYCLES  = 

I  = 

XU)  = 

YU)  = 

4000. 

1 

0 

172 

2 

062 

141 

.085 

.116 

4 

.0SE 

.056 

5 

.108 

.073 

6 

.116 

.051 

7 

.126 

.028 

a 

.130 

.016 

9 

.128 

0 

CYCLES  = 

I  = 

X(I)  = 
Y(I)  = 

5000. 

1 

0 

.178 

2 

.039 

.163 

3 

.064 

.147 

4 

.084 

.130 

5 

.109 

.034 

6 

.119 

.063 

7 

.125 

.039 

8 

.129 

.011 

9 

.128 

0 

CYCLES  = 

I  = 
X(I)  = 
YU)  = 

6000. 

1 

0 

.133 

2 

.031 

.188 

3 

.055 

.178 

4 

.083 

.154 

5 

.103 

.126 

6 

.120 

.097 

7 

.128 

.058 

8 

.134 

.015 

9 

.129 

0 

CYCLES  = 

1  = 
XU)  = 
Y(I)  = 

7000. 

1 

0 

.1SS 

2 

.061 

.183 

3 

.088 
.  152 

4 

.115 

.118 

5 

.125 

.093 

6 

.131 

.076 

7 

.134 

.051 

3 

.135 

.034 

9 

.136 

0 

CYCLES  = 

I  = 

XU)  = 
YU)  = 

9000. 

1 

0 

.213 

2 

.039 

.203 

.060 

.200 

4 

.087 

.175 

5 

.115 
.  146 

6 

.125 

.127 

7 

.137 
.  0S5 

8 

.141 

.061 

9 

.139 

0 

CYCLES  = 

I  = 
XU)  = 
Y(I)  = 

10000. 

1 

0 

.224 

2 

.035 

.218 

.064 

.204 

4 

.083 

.183 

5 

.113 

.160 

6 

.125 
.  133 

7 

.136 

.101 

3 

.142 

.065 

9 

.137 

0 

CYCLES  = 
I  = 

XU)  = 
Y(I)  = 

110C0 

1 

0 

.238 

2 

.  032 
.229 

3 

.060 

.217 

4 

.083 

.158 

5 

.113 

.167 

6 

.  135 
.125 

7 

.141 

.093 

3 

.143 

.065 

9 

.137 

0 

CYCLES  = 
I  = 

X(I)  = 
YU)  = 

12000 

i 

.001 

.242 

2 

.042 

.233 

.063 

.223 

4 

.091 

.205 

S 

.113 

.180 

6 

.131 

.155 

7 

.143 

.121 

3 

.151 

.067 

9 

.145 

0 

CYCLES  = 

I~  = 

:<U)  = 

Y(  I)  a 

13000 

C 

.245 

.CIS 

.251 

•J 

.cs: 

,22  5 

4 

.CSC 

.223 

Use 

S 

.  140 

.158 

7 

.153 

.116 

3 

.131 

.063 

9 

.153 

0 

118 


Table  8  continued 


CYCLES  =  14000. 

I  ®  1 

xd)  ■  o 

Yd)  =  .270 

2 

041 

264 

3 

.077 

.245 

4 

.105 

.228 

5 

.134 
.  1S6 

6 

.150 

-.163 

7 

.161 

.129 

8 

.168 

.079 

9 

.161 

0 

CYCLES  =  15000. 

I  *  1 

Xd)  a  0 

YCI)  a  .283 

2 

044 

276 

3 

.079 

.261 

4 

.112 

.234 

5 

.139 

.203 

6 

.160 

.162 

7 

.173 

.117 

8 

.177 

.070 

9 

.163 

0 

CYCLES  =  1S000. 

1  a  1 

Xd)  =  0 

Yd)  =  .293 

2 

,049 

.285 

3 

.091 

.264 

4 

.120 

.239 

5 

.146 

.205 

6 

.161 

.177 

7 

.177 

.138 

8 

.133 

.088 

9 

.167 

0 

CYCLES  =  17000. 

I  *  1 

Xd)  *  0 

Y(I)  =  .319 

2 

.042 

.314 

3 

.OSS 

.291 

4 

.120 

.255 

5 

.162 

.211 

6 

.177 

.180 

7 

.183 

.141 

a 

.196 

.085 

9 

.179 

0 

CYCLES  =  18000. 

I  =  1 

XCI)  a  0 

Yd)  =  .338 

2 

.062 

.333 

3 

.105 

.303 

4 

.133 

.274 

5 

.162 

.233 

6 

.181 

.193 

7 

.196 

.157 

3 

.205 

.084 

9 

.186 

0 

CYCLES  =  19000. 

I  =  1 
X(I)  a  0 
Y(I)  a  .351 

2 

.065 

.350 

3 

.103 

.323 

4 

.127 
.  295 

5 

.157 

.262 

6 

.187 

.208 

7 

.203 

.167 

3 

.214 

.108 

9 

.204 

0 

CYCLES  a  20000. 

I  =  1 

Xd)  a  0 

Y(I)  =  .386 

2 

,058 

.376 

3 

.111 

.349 

4 

.156 

.302 

5 

.186 

.251 

6 

.207 

.207 

7 

.221 

.169 

3 

.232 

.092 

9 

.221 

0 

CYCLES  =  21000 

I  =  1 

Xd)  =  0 

YCI)  =  .408 

2 

.061 
.  395 

3 

.113 

.356 

4 

.  158 
.321 

5 

.201 

.251 

6 

.225 

.190 

7 

.239 

.142 

3 

.245 

.081 

9 

.232 

0 

CYCLES  =  22000 

I  a  1 

Xd)  =  0 

Y(I)  =  .444 

2 

.080 
.  432 

3 

.  135 
.387 

4 

.176 

.342 

5 

.213 

.281 

6 

.237 

.228 

7 

.252 

.167 

3 

.253 

.039 

9 

.249 

0 

CYCLES  a  =25 CC 
I  =  1 

X(I)  =  C 

Yd)  =  .475 

i. 

2 

.  C7S 
.<•71 

.123 

.431 

4 

.170 

.384 

5 

.208 
.  333 

6 

.241 

.264 

7 

.262 
.  139 

3 

.274 

.105 

9 

.262 

0 

119 


Table  8  continued 


CYCLES  = 

I  = 

23000. 

1 

2 

3 

4 

5 

5 

7 

3 

9 

XC 15  = 

0 

.  0G4 

.113 

.151 

.208 

.245 

.259 

.284 

.26? 

Yd )  = 

.4S5 

.453 

.45S 

.414 

.349 

.280 

.204 

.112 

0 

CYCLES  = 

I  = 

23500. 

1 

2 

3 

4 

5 

6 

7 

3 

9 

XC  1 5  = 

0 

.077 

.145 

.201 

.252 

.274 

.285 

.293 

.279 

YC 1 5  = 

.524 

.50? 

.457 

.  389 

.307 

.244 

.195 

.126 

0 

CYCLES  = 

T  = 

24000. 

1 

2 

3 

4 

5 

6 

7 

3 

9 

XC  1 5  = 

0 

.081 

.139 

.208 

.249 

.272 

.295 

.308 

.295 

YC  1 5  = 

.552 

.540 

.498 

.416 

.341 

.  cSO 

.216 

.117 

0 

CYCLES  = 

I  = 

24300. 

1 

2 

3 

4 

5 

6 

7 

3 

9 

XC  15  = 

0 

.087 

.134 

.128 

.241 

.279 

.306 

.314 

.299 

YC  15  = 

.575 

.553 

.520 

.464 

.383 

.300 

.203 

.131 

0 

CYCLES  = 

I  = 

24500. 

1 

2 

3 

4 

5 

6 

7 

3 

9 

XC  1 5  = 

0 

.09? 

.149 

.  2C6 

.231 

.285 

.312 

.313 

.304 

YC  1 5  = 

.579 

.567 

.522 

.447 

.373 

.299 

.205 

.119 

0 

CYCLES  * 

I  = 

24750. 

1 

2 

3 

4 

5 

6 

7 

3 

9 

XC  1 5  = 

.001 

.091 

.149 

.190 

.251 

.290 

.313 

.327 

.314 

Y(I)  = 

.600 

.580 

.534 

.491 

.362 

.309 

.224 

.115 

0 

CYCLES  = 

I  = 

25000. 

1 

2 

3 

4 

5 

6 

7 

3 

9 

XC  1 5  = 

0 

.063 

.133 

.191 

.237 

.276 

.310 

.  336 

.319 

Y(I)  = 

.628 

.615 

.575 

.511 

.44? 

.373 

.282 

.125 

0 

CYCLES  = 

I  = 

25250. 

1 

2 

3 

4 

5 

6 

7 

3 

9 

Xd)  = 

0 

.CSS 

.149 

.230 

.  284 

.319 

.334 

.342 

.318 

Y(I)  = 

.648 

.628 

.583 

.470 

.374 

.  280 

.190 

.101 

0 

CYCLES  = 

I  = 

25500. 

1 

2 

3 

4 

5 

6 

7 

3 

9 

Xd)  = 

0 

.085 

.154 

.219 

.283 

.295 

.324 

.351 

.339 

Yd)  = 

.665 

.648 

.604 

.523 

.447 

.372 

.289 

.126 

0 

CYCLES  = 

25750. 

1 

2 

5 

3 

7 

3 

9 

Xd)  = 

C 

.CSC 

.153 

•  sc? 

.255 

•  oC  0 

.  335 

.332 

.347 

Yd)  = 

.687 

.  See 

.  638 

.  S7i 

.  -S3 

.  4G0 

.304 

.122 

0 

CYCLES  = 

I  = 

26000. 

1 

2 

3 

4 

C 

6 

7 

8 

9 

XC  I )  = 

0 

.082 

.ISO 

.  23S 

•  cSc 

.338 

.359 

.375 

.353 

Yd)  = 

.703 

•  (  0  3 

•  636 

.577 

.4S5 

.358 

.273 

.129 

0 

120 


22260- 
32/  S  0  - 

32;  OC  - 
32050  — 
31000  — 
3l8oo  — 
3/800  - 

31  700 — 
3/  600 

3)800 


I 


Table  9  Digitized  Data  for  Snow  Test  9 


NUtt  a  12  SR  3  .743  W  a  7.970 
T  a  .704  LOAD  =  3400. 


CYCLES  = 

I  = 

31300. 

1 

2 

3 

4 

5 

6 

7 

3 

9 

X(  I  5  a 

0 

.112 

.168 

.215 

.260 

.250 

.323 

.366 

.352 

Y(  I  )  = 

.S7S 

.647 

.599 

.540 

.478 

.415 

.335. 

.117 

0 

CYCLES  a 

I  = 

31400. 

1  2 

3 

4 

5 

6 

7 

8 

9 

XCI)  = 

0 

.085 

.157 

.217 

.272 

.315 

.354 

.373 

.361 

Y(  I  )  = 

.639 

.672 

.613 

.545 

.455 

.374 

.246 

.120 

0 

CYCLES  a 

I  = 

31500. 

1  2 

3 

4 

5 

6 

7 

8 

9 

XCI)  = 

0 

.069 

.135 

.213 

.274 

.315 

.353 

.375 

.365 

Y(  I  )  = 

.692 

.678 

.643 

.553 

.462 

.383 

.271 

.111 

0 

CYCLES  a 

I  = 

31600. 

1  2 

3 

4 

5 

6 

7 

8 

9 

XCI)  = 

0 

.089 

.179 

.233 

.277 

.310 

.355 

.390 

.377 

VC  I  J  = 

.705 

.651 

.618 

.S54 

.484 

.419 

.316 

.107 

0 

CYCLES  a 

I  = 

31700. 

1  2 

3 

4 

5 

6 

7 

8 

9 

XCI)  a 

0 

.114 

.168 

.237 

.306 

.350 

.375 

.390 

.375 

Y(  I  )  = 

.716 

.681 

.634 

.547 

.434 

.339 

.232. 

.096 

0 

CYCLEj  a 

31800. 

1  a 

1  2 

3 

4 

5 

6 

7 

8 

9 

XCI)  a 

,002  .070 

.144 

.218 

.268 

.312 

.366 

.389 

.376 

V  <  I  3  a 

.730  .712 

.667 

.584 

.512 

.433 

.270 

.117 

0 

CYCLES  a 

I  a 

31900. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

XCI)  * 

0 

.124 

.181 

.249 

.310 

.346 

.378 

.395 

.378 

YC  I  )  a 

.725 

.631 

.634 

.545 

.441 

.360 

.231 

.111 

0 

CYCLES  » 

I  « 

32000. 

1 

2 

3 

4 

5 

6 

7 

3 

9 

XCI)  a 

0 

.079 

.185 

.  24S 

.309 

.346 

.381 

.403 

.389 

YC  I )  » 

.738 

.722 

.640 

.564 

.459 

.380 

.282 

.124 

0 

CYCLES  » 

I  ■ 

32050. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

XCI)  a 

0 

.080 

.183 

.251 

.307 

.359 

.388 

.407 

.335 

Y  (  I  )  a 

.746 

.740 

.663 

.580 

.482 

.376 

.266 

.116 

.014 

CYCLES  a 

I  a 

32100. 

1  2 

4 

5 

c 

7 

8 

9 

XCI)  a 

0 

.  120 

.  1  £4 

.256 

.307 

.  3<S 

.330 

.408 

.395 

YC  I )  a 

.748 

.710 

.654 

.563 

.477 

.  3S7 

.300 

.127 

0 

122 


Table  9  continued 


CYCLES  *  32150. 

1=123456789 
X(I)  ■  0  .099  .167  .228  .289  .344  .376  .415  .400 
Y(I)  -  .761  .729  .683  .627  .528  .424  .336.  .128  0 


CYCLES  =  32200. 

1*123456789 
XCI)  *  0  .128  .205  .265  .313  .351  .388  .418  .405 

VC  1 3  *  .770  .  724  .  647  .  573  .  495  .  415  .  296  .126  0 


123 
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